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Abstract 

The low-energy properties of the two-dimensional Heisenberg model with spin-| on 
a square lattice are investigated on the basis of the local dimer order. The lattice is 
divided into square blocks consisting of the quartet of spins. The spin variables and the 
Heisenberg Hamiltonian are expressed in terms of the low-energy quartet variables. On 
the basis of the Dyson-Maleev representation the spin-wave theory of the quartet state 
is developed. The spectrum of the lower magnon excitations consists of three degenerate 
modes with the energy gap A = 0.17 J. The ground state energy per spin E/N = —0.6 J. 
This preprint repeats in the main the previous one but it contains calculations of the basic 
corrections and therefore has complete character. 
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I. Introduction 



Unusual states of the two-dimensional Heisenberg antiferromagnet attract attention in con- 
nection with the problem of the magnetic state of the cuprate superconductors [p]] and very 
complete references in this review. This activity was stimulated by the experiments which 
demonstrated that superconductivity in cuprates is realized in the strongly correlated param- 
agnetic spin state which was named the spin-liquid state. 

Many different models of the spin-liquid state were proposed [jj - The spin-liquid state 
as a linear combination of different dimer states was considered in the works || |J and was 
named the RVB-state. The dimer represents the state of the two spins-| with the total spin 
equal to zero. The dimer can be formed from two neighboring spins- 1 as well as from two 
separated spins- ~. Numerical simulation of the RVB-state was produced in the work || and 
it was demonstrated that this state for the two-dimensional Heisenberg antiferromagnet with 
the nearest neighbor exchange had very low energy. This energy is higher but very close to the 
energy of the Neel state. Of course, the main problem is the transformation of the magnetic 
state with doping. But for progress in this problem we must have sufficiently simple model of 
the spin-liquid state because the original RVB-state permits only numerical consideration. 

Such simple model of the spin-liquid state was proposed in the works || [E(J. The plane 
is split into square blocks containing four spins. The block size is 2a, where a is the lattice 
constant. The complete set of eigenstates for such block or quartet consists of 16 states which 
can be easily found for the exchange interaction J\ between the first and J 2 between the second 
neighbors. For J 2 < J\ the ground state \ip > of such quartet is the state with the total 
momentum equal to zero. This state can be presented as a sum of the two possible dimer states 
for the nearest spins. Its energy is — 2J X + J2/2. 



As it was proposed in [0, [L(| that one can get good approximation for the spin-liquid state 



if \ip > state is chosen as ground state in each quartet and the quantum fluctuations are taken 
into account. At J 2 = the energy per site was in [[UJ E = — 0.655J that is very close to the 
numerical results |4j] for Eq = — 0.668J. In other paper || based on the quartet approach some 
sort of the Schwinger boson representation for the Hubbard operators of a quartet was used. 
In this work the energy was sufficiently high E Q = —0.57 J. The correspondence between the 
numerical simulation and the quartet model is not so obvious, because in the quartet approach 
the quantum fluctuations are taken into account. But approximations of || |TU[ were sufficiently 
crude. The contribution of the higher-energy quartet states to the energy of the ground state 
and to the spin dynamics was omitted. In the work [10] uncontrollable approximation based 
on some decoupling scheme for the Green's functions was used. In the work || the Schwinger 
boson method is applicable only over parameter 1/N and it is also sufficiently crude. 

In this work, which is based on the main idea of [|| [TO] about the quartet ground state for the 
Heisenberg antiferromagnet, we use spin- wave approach to description of the triplet excitations. 
All approximations of this work are controllable. In this work we restrict our consideration to 
the case of the unfrustrated Heisenberg antiferromagnet with J 2 = 0. We consider contributions 
of the lower singlet state and all triplet states into statistics and dynamics of the quartet state 
of the Heisenberg antiferromagnet. The role of the additional singlet state and the quintet state 
is essentially less because direct transitions from the lower singlet state into these states are 
absent. We do not consider these states and it will be evident that their contribution into the 
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ground state energy is sufficiently less. 



II. The effective Hamiltonian of the model 

The Hamiltonian of the Heisenberg antiferromagnet has a well-known form: 

H=J J £(*i*i>) (2- 1 ) 

<i,v> 

where < IV > notes summation over the nearest neighbors. For construction of the quartet 
magnetic state let us divide the square lattice into four-spins blocks ||, f[(| 

< 2 - 2) 

and find all eigenstates of the quartet of spins (see Appendix A). Every quartet has 16 states: 
two singlets \<f >, \ip > , three triplets \tm >, \xm >, \ym > and quintet \qm >. If we restrict 
our consideration to the low-energy singlet ip and triplet t then the spins S; consisting a quartet 
can be expressed in terms of the Hubbard operators acting in \ip >, \tm >, m = ±1, subspace 
at the quartet n: 

= \<pn >< rup\, Z* m '* m ' = \tmn >< nm't\, 
Z^ tm = \ipn >< nmt\, = \tmn >< rup\, (2.3) 

The spin operator s ni consisting of the quartet with the number n at the doubling square lattice 
may be represented as a sum of two vectors: 

s ni — ( — )* +1 Lj„, + -Sun (2-4) 

where L tn mixes \<p > and \tm > states 

L ; n = ^=(Z^° + Z^), L^ = 1={Z^- X -Z^), L m = ^=(Z^-Z^) (2.5) 
for p = t and S u , n are operators of the spin 1 acting in t-subspace 

dpqn — ^n J ^pqn ~ V 6\A n + ^ n ), O pqn — V A\& n + ^ n ), 

for p,q — t, where A ± = A x ±iA y are the spherical components of a vector A. Using expression 
of (|A.3|) in Appendix A for the energies of \(p > one can get the Hamiltonian of the problem 



tin 



-(J/3) Yl (2Z n m '*Z% tm + (-y^z^z^ + zf m z^- m ) 

<nn'>m 

+ (J/8) Z n n ' tm ' (zT'' tm - (-) m+m ' Z^-" 1 ') . (2.7) 

<nn'>mm' 
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The Hubbard operators Z^ b for a,b = ip, tm are not convenient for the further analysis of the 
Hamiltonian ( |2.7| ). Therefore we map the Gilbert space consisting of \(p > and \tm > states 
into the Gilbert space of the three-dimensional Heisenberg algebra t^,t m . Where t+(£ m ) are 
creation (annihilation) Bose operators. The general state of such algebra is |n + ,n_,no > and 
mapping has simple form 

|v> |o,o,o >, |*i >-► |i,o,o>, 

|t-l >-► |0,1,0>, [*0 > — »• |0, 0, 1> (2.8) 



This mapping generates the Holstein-Primakoff representation |TTJ for our Hubbard operators 



(z^)hp = p(i - N n )p, (z% tm ) HP = pyi - N n t nm P, 

(Z^ ,<p )hp = Pt\ m \j\ - N n P, (Z* m '* m ) HP = Pt+ m t nm rP, 

N n = J^nm^nm; (2-9) 



where P is a projector on the lower four states (|2.8|) of the Heisenberg algebra £+ m , t nm . At 
the low temperatures T <C J or a small excitation of oscillators \nm > one can produce an 
isometric transformation at each site 

{Z°" b ) DM = V-\Z a ' b ) HP V (2.10) 

and omit the projector P. One can choose the operator V from the condition of elimination of 
roots in ( [2. 10 ), and get the Dyson- Maleev representation for the Hubbard operators [12, |H| 



(Z^)dm = (1 - N n ), (Zf m ) DM = t r , 

{Z™ ,(P )dM = tnmO- - N n ), (Z* m >* m ) DM = 

N n — '^^t nm t nm . (2- 11, 



nmi 
+ 4- 

nm L nm' i 



The Dyson-Maleev representation for Z-operators is not a direct operator identity but it can 
be used at low temperatures and a small level of excitations. Substituting the Dyson-Maleev 
representation ( |2.11| ) in the Hamiltonian ( |2.7| ) we get the following effective Hamiltonian 



H = -(NJ/2) + JJ2N n -(J/3) J2 (2t+ m (l - N n )t n , m 

n <nn'>m 



+ (") m (4m(l - NnK_ m (l ~ N n >) + Wn'-m) ) 

+ (^/ 8 ) Yl ttrJnm' (*n'm'*n'm ~ (-) m+m t+,_ m t n >- m ') , (2.12) 



<nn'>mm' 



where N is the number of spins in the plane. This Hamiltonian is nonhermitian due to properties 
of the Dyson-Maleev representation ( |2.11| ) but if we use it for the computation of Green's 
functions any contradictions do not appear. 
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III. The approximate solution of the effective Hamil- 
tonian 



If we consider only the quadratic terms in the Hamiltonian Q2.12Q we can easily diagonalize 
them with the help of u — v transformation (these quadratic terms are Hermitian) and get the 
square of the magnon energy in the momentum representation 

El = J 2 (l-8 7k /3), 7 k = (cos^a) + cos(2A;,a))/2, (3.1) 

where a is the distance between spins, k is the quasimomentum. One can easily see that E\ is 
negative at 7k > 3/8. This means that our system is unstable in the quadratic approximation 
and the higher terms in the Hamiltonian ( 2.12 ) must be taken into consideration for solving the 
problem of stability. We will treat the Hamiltonian ( |2.12| ) by the mean-field method. The same 
can be done by the Green's function method. The results are identical. We seek the solution 
which is invariant over the rotation and consider the following normal and abnormal averages: 

< tkmtk'm' >= 5 m ,m"5(k - k')iV k , 

< ttji'm* >= R m 5 m ,- m '5(k + k')F k +, 

< ^km^k'm' >= 5rn,- m > 5 (k + k )F k . (3.2) 

Our approach based on the Dyson- Maleev representation ( [2.11| ) is correct if the contribution of 
the higher states at every lattice sites is small. Because of that the following conditions must 
be fulfilled 

{N p ,Ff) = E(7k) P (iV k ,F k ± ) « 1, (3.3) 

k 

where the sum over k notes the normalized integral over the Brillouin zone. In this approx- 
imation we can omit the contribution of the members of the six order and get the following 
mean-field Hamiltonian: 



H m f - E mf + (ak^km^km + (-^(Ac^km^-k-m + ^km^-k-m) /2J , (3.4) 



km 



where E m j is the initial ground state energy in the mean-field approximation: 

E° mf = (N J/2) ( - 1 + (3# - m ~ 2Ff - 2F+)/4 

-(3/8) ]T(2a k iVk + faF+ + 7^)) (3.5) 

k 

and the coefficients ak, /?k, <5k are 

a k = 4(3/4 + 4(JVi + F+) + (4iV + F+ - 1 + 3iY 1 /4) 7k )/3, 
A< = 4(2(JVx + Ff) + (8iV - 1 - 3Ff /4) 7k )/3, 

4 = 4(2F+ - 1 - 3F+/4) 7k /3. (3.6) 

The quadratic part of the Hamiltonian H m f ( p.4|) was obtained from H fl2.12p by computing 
all possible paired averages in the fourth order terms of H (|2.12 ) . The initial energy E m ^ is 
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determined from the condition of coincidence of the average values of the Hamiltonian (p. 12 ) and 



This procedure corresponds to determination of the one-particle energy by the Green's 
function method in the one-loop approximation. The Hamiltonian ( |3.4|) can be diagonalized 
with the help of the u — v — w transformation (this procedure is equivalent of solving the Dyson 
equation for normal and abnormal Green's functions): 



tkm — "k&km + (~ ) m,y k&_ k _ m! 
+ 



t: km = U kKm + (-) m W k 6_ k _ m , 

«k - ^k^k = 1- (3.7) 



After this transformation the Hamiltonian (13.41) takes the form 



H m f = E m f + ^ Ekb^ m bk m (3.8) 

km 

with the ground state energy 

E mf = E° mf + (3/4)iy£(£ k - J Qk ) (3.9) 

k 

and 



E k = JR k , R k = ^al - /3 k 5 k . (3.10) 
The coefficients ii k , v k , w k have form: 



«k = \ ("k + Rk)/2R k , f k = z k /3j 



k- 



Zk = — \f\ock - -Rk|/2|/3 k 5 k |-Rk, w k = 2 k 5 k . (3.11) 

So far as the coefficients a k , /3 k , # k depend on the overages N p , and consequently the coeffi- 
cients w k , v k , u> k also depend on these parameters we can get the closed system of equations if 
we substitute Eqs. ( |3.7|) for t km , t^ m into Eqs. (|3.2| ), ( |3.3|) which determine N p ,F^. 

iVk= K/i? k )K+l/2)-l/2, 

F± = -[(5 k ,/? k )/i? k ](n k + 1/2), (3.12) 

where n k is the Planck function: n k = (exp (j3Ek) — l)~ x and /3 = 1/T is the inverse temperature. 
The system of equations of Q3.12 ) was solved numerically (see Appendix B for details) and we 



get the averages N p , as a function of the temperature. We also computed the free energy 
of the quartet state from ( |3.8|) and determined the magnon energy i? k 



£ k = VA + 2Be + Ce 2 , e = 7k 
£ k w A(l + k 2 £ 2 ), at ka < 1 

A = VA + 2B + C, £ 2 = -(£ + 2C)a 2 /8A 2 (3.13) 
All these parameters are presented in Table 1. 
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T/J 


F/NJ 


A/J 


A/J 2 


B/J 2 


C/J 2 





-0.567 


0.174 


3.474 


-1.595 


-0.254 


0.1 


-0.567 


0.205 


3.484 


-1.592 


-0.258 


0.3 


-0.570 


0.416 


3.525 


-1.562 


-0.229 


0.5 


-0.585 


0.704 


3.452 


-1.412 


-0.133 



T/J 


N 


iVi 


F + 


Ft 




Fr 





0.045 


0.025 


0.087 


0.137 


0.035 


0.083 


0.1 


0.046 


0.026 


0.087 


0.137 


0.034 


0.083 


0.3 


0.050 


0.029 


0.086 


0.136 


0.028 


0.076 


0.5 


0.068 


0.030 


0.079 


0.131 


0.011 


0.055 



TABLE 1. The result of calculation of the free energy F, the energy gap A, the parameters of 
the magnon energy A, B, C and the normal and abnormal averages N p , F^. 

We can see that all normal and abnormal averages are small and do not exceed 0.14. This 
property of our solution justifies dropping of the six-order term in the effective Hamiltonian 
(|2.12j) . The smallness of these averages justifies also using of the Dyson- Maleev representation 
for the Hubbard operators. Really, one can check that the occupancy W n+n _ no of the state 
|n + , n_, n > of the Heisenberg algebra t+, t m for our u — v — w transformation ( |3.7|) is 
determined by the formula 

W n+n _ no = (1 + iV )- 3 (j^-J (3-14) 

and the total occupancy of the unphysical states W unp is equal to 

W unp w (6iV 2 + 4iV 3 + iV 4 )/(l + iV ) 4 = 0.044. (3.15) 
All this prove that we found correctly the rotationally invariant ground state of the Hamiltonian 

IV. Corrections to the mean-field energy of the ground 
state. 

For checking accuracy of our solution we compute the second order correction to the ground 
state energy in the mean-field approximation ( |3.9| ). This correction is presented by the Feynman 
diagram in Fig.l. 
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Fig.l The Feynman diagram for the second order correction to the energy of the ground 
state. 

We restrict our consideration to the case of the zero temperature. The part of the Hamilto- 
nian ( [2.12j ) essential for computation of the correction of Fig.l after u-v-w transformation has 
form 

Hm = (^/12) {[IGUkiUka^kaTkifak* + ^k 4 ) + 

ki,m,n 

37k 1+ k 2 M kl f; k2 (M k3 f; k4 - u^v^-^^b+^b+^bt^bt^J + 

(frmk ~+ bmk, U k -> W k , V k -> Uu)}. (4.1) 

Here summation ki, k2, k3, k4 is produced over the surface ki + k2 + k3 + k4 = 0. 

The first term in the square brackets in (|4.1|) follows from the L — L interaction and the 
second term follows from the S — S interaction. The last term in the figured brackets is the 
conjugated to the first term. On the base of the Hamiltonian one can easily compute the second 
order correction of Fig.l to the ground state energy. We compute the six-dimensional integral 
over ki, k2, k3 by the Monte-Carlo method and getting the following result SE = 0.005NJ. 

This correction to the ground state energy is sufficiently small that verifies correctness of 
our mean-field solution of the Hamiltonian ( 2.12| ). The validity of the perturbation theory is 



based on the smallness of angle integrals over k because another numerical small parameters 
are absent in our theory. In fact, the small parameter is 1/z where z is the number of the 
neighboring spin to our quartet of spins. 

V. Influence of higher energy levels of the quartet on 
the ground state energy 

It is clear that transitions at the higher energy levels of the quartet reduce the ground state 
energy of the quartet state of the two-dimensional Heisenberg antiferromagnets. As was noticed 
above the most essential contributions originated from the higher triplet states \xm >, \ym > 
(see (|A.5|) . In order to consider this additional to \tp > and \tm > states we must modify 
Eqs. (|2.4j )- (|2.6| ) for the spin operators Sj taking into account contribution of \xm >, \ym > 
states. On the base of the quartet wave function ( A.5|) we get the following representation for 
the spin operators Sj 



4s 4 = 4(-)* +i L t - 2(L X - (-yL y ) + ^ + + S yy 

+2t(Ste + S xt — (— Y(Sty + Syi)) + (— Y(s xy + S yx ) (5-i, 
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where z\ = (1,1,-1,-1) for i = 1,2,3,4. The operator of the triplet excitations L p and 
the spin-1 operator S pq for p,q — t,x,y; m = 0, ±1 are expressed by the same manner as 
in Eqs.( |2.5|) , ( |2.6|) . For the further analysis we replace the Habburd projection operators 
Z^, Z%' pm , ZP" 1 ^, ZP m ' qm ' by the more convenient Bose operators x+ m , y+ m , t+ m {x nm , y nm , t nm ) 
creating (annihilating) bosons of x, y, t type with projection m at site n. The corresponding 
Dyson- Maleev representation for the Hubbard operators is of the form ( 2.11| ). 

(zr) D M = (i-N n ), ^ (z^ m ) DM = Pnm , 

(ZP m ' v ) DM = p nm (l - N n ), (^z pm ' qm ) DM = p nm q nm ', 

Nn = Y,PtmPn m - (5-2) 



Hamiltonian ( |2.1|) can be represented as 

H — J y^(si n + S3 n ) (S2n + S4n) + J /j (^2n^3n' x + ^ln^4n' x ) + 

n <n,n' x > 

J 2l ( S 2n s ln^ + s 3n s 4n.^) (5-3) 

<n,n' y > 

where < n,n' x > (< n, n y >) means summation over the nearest neighbor of block n in x (y) 
direction. Using the representations ( |5.1| ),( pT2"D for the spin operators s, one can rewrite the 
Hamiltonian ( |5.3|) in the momentum representation: 

H = H t + H x + H y + H tx + H ty , 

H x = JJ2(( 2 + ^k/3)Xk m x km + (-) m rik(x£ m x\_ m + x km x_ k _ m )/6) , 

km 

H V = J E (( 2 ~ W%kmJten» - (-) m 1k(«/km!/-k- m + VUmV -U-m) / ^) , 
km 

H tx = (zJ/3) J2(-) m -M2k x a) (t km + (-) m t\_ rn ) (x^ m + (-) m x+ m ) , 

km 

H ty = -(iJ/3) £(-) m «n(2V) (t km + (-) m t\_ m ) (y^ m + (-) m yL) • (5-4) 



km 



where r/ k = (cos(2k x a) — cos(2k y a))/2 and H t is the Hamiltonian (|2.12j ), which was discussed 
previously. In H x , H y , H tx , H ty we excluded from our consideration all terms more than the 
second order by creation and annihilation operators. These terms lead in the mean-field ap- 
proximation to some renormalization of coefficients in the quadratic form (|5.4j ) which does not 
effect essentially on results obtained below. 

The further plan is as follows. We take into consideration the Hamiltonians H x and H y 
precisely and consider the Hamiltonians H tx and H ty as perturbation. We will check that such 
approach is suitable for the computation of corrections to the ground state energy. 

The Hamiltonians H x and H y can be easily diagonalized with help of the u—v transformation 
The coefficients a k ,/5 k and 5 k in this case are 

a k = J (2 ± 7? k /3), p k = 5 k = ±J?? k /3, (5-5) 
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where the upper sign "+" corresponds to H x and the lower sign "— " H y . As a result we get 
the energy for the x, y quasiparticles in a form 



^ = \K - PI = 2 J^Jl ± W3 ~ J(2 ± ry k /3 - r^/36). (5.6) 
The average number of quasiparticles x, y in the ground state is sufficiently small 

< n x >=< n v >= VdkR; — V ril = — < N . (5.7) 
y Y 144 Y 576 V ; 

Hence, the contribution of H x and if,, in the ground state energy can be estimated as 1CT 3 J 
per spin. 

Notice, that dependence of the energies e^ v on k as it follows from Eq. (|5.7|) is sufficiently 
weak. This leads to the strong interaction inside the x — y sector of the quartet Gilbert space. 
As it follows from Eq. (|5.1|) there are nonlinear terms in the total Hamiltonian of the spin-spin 
interaction which lead to the strong interaction inside of the x — y sector. These terms in the 
Hamiltonians have a characteristic form: 

H mt = (2J/V3)J2 7k(Lt i k(S Z j /j -k + Syx-k) (5.8) 

k 

This means that, in fact, we do not know the structure of excitations in the x — y sector. 
But for computation of corrections to the ground state energy we can use the zero-dimensional 
approximation for the Green's function of the x, y quasiparticles fl4| . In this approximation 



the energy of x and y particles is simply 2J ( |5.6|) and their Green's function can be easily 
found. One can show that in result of the interaction ( |5.8j ) a form of the Green's function in 
the energetic space became the Gauss form instead the Lorentz form as in the case of the weak 
interaction. But what is the most important for us - a position of the center of peak in the 
energetic space of the imaginary part of the Green's function does not change at e = 2 J. The 
corrections to the ground state energy due to the Hamiltonian ( |5.8| ) are equal to zero in this 
approximation, because a contribution of the condensate the x — y particles are small Eq. (|5.7|) . 

The Hamiltonians H tx and H ty lead to an equal contributions to the ground state energy, 
which can be calculated in the second order of the perturbation theory: 



dE mf =-(A//2)2^ ■ 

qZ fc q 

= -( J 2 iV/6) E *m 2 (2M K + + = -0-024JiV, (5.9) 



£q + £q 



where u, v, w are defined in ( 3.1 1|) and iV is the number of spins in a plane. We see that the 



interaction of x and y particles with t particle decreases the ground state energy by 0.024J per 
spin. This correction will be reduced if we take into consideration excluded from Hamiltonian 



( |5.4|) terms of third and forth order over the creation and annihilation operators. 
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VI. Properties of the solution and discussion 



The energy per spin of the quartet state at the temperature equal to zero is sufficiently high: 
E/N = — 0.6J. The coefficients B and C in the magnon energy (|3.14|) are negative and the 
energy has minimum at k = and possesses a gap. This gap is equal to 0.174 J at T = 0. 
Because our excitations have a gap, any long-range order is absent in the quartet magnetic 
state. We can easily calculate the simultaneous spin correlation function using representation 
for spin operators and the representation ( 2 . 1 1|) for the Hubbard operators: 



Kt'n'i' =< S ni S n'i' > = <W-Kii'( r ): ^ = T n - T n , , (6.1) 

where i, i = 1, 2, 3, 4 numerate spins inside quartet, a, a' = 1,2,3 are vector indices. According 
to ( j2.4|) we have 

K iV {v) = (-)*+*' < n n n n , > +1 < sffi > . (6.2) 



Substituting ( |2.5| ), ( |2.6| ) in ( |6.2j ) and using ( |2.11| ) one can get 



K w {r) = (2(1 - F + - 4iV )iV(r) + (1 - 2F +)F-(r) + (1 - 8iV )F + (r) 

+l(iV 2 (r)-F-(r)F + (r)), (6.3) 

were N(r), F ± (r) are the Fourier images of correlators N^, F^: 

(N(r),F±(r)) = 5> X p(-ikr)(iV k , (6.4) 



When we derived (|6.3| ) we take into consideration the second and fourth order over t , t terms. 
The contribution of the six-order terms is small over N p , Ff. The Ion g-range order is absent in 
our quartet state and asymptotic behavior is determined by the gap in the magnon spectrum 
A and the effective magnon mass. It can be easily found from ( |6.3|) , ( |3.12|) , ( p^6|) and ( |3.14|) : 



Kij{r) = J D(-) i+J (a/0 2 (27rrO" 1/2 exp(-r/0, r>^ (6.5) 

The constant D is a function N p , F^ and the correlation length £ is determined by the magnon 
energy. Taking into account interaction of x and y particles with t particle does not lead to 
the efficient modification of the correlation function since this interaction is strongly reduced at 
small k due to the structure of contributions of the x, y particles into the correlation function 
which are proportional to sin(2k x a) and sin(2k y a). Therefore, a corrections to Kij is of the 
order of SK^ ~ (a 2 /C, 2 )Kij <C K^. The behavior of the correlation function at the large 
distance is similar to an antiferromagnet with short-range order |TJ, |7|, [5] and contains the 
"stagger" factor (— . 

Our results essentially differ from results obtained in works H, [TO] where the quartet state 



was proposed. We considered only the two-dimensional unfrustrated Heisenberg antiferromag- 
nets. Our consideration is quantitative and does not contain unjustified assumptions. We did 
not find any gapless excitations as in || [|, |10[ . Moreover our magnetic excitations have the 
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IS 



similar structure as in |Hj, but the energy is quite different. The ground state energy in flC 
sufficiently low E /N = -0.655J. 

We connect these contradictions with the crude method of solution of the Hamiltonian 
(f2.7|). The method of decoupling of Green's functions used in [[10] is rather indefinite and result 
depends on the method of decoupling. It is reasonable to suppose that low energy obtained 
10] as well as in the Schwinger boson method |J is an accidental result of approximation. 



in 



The corrections to these approximations are not small and can destroy this low energy. We 
believe that application of the Schwinger boson method to the Hamiltonian ( p.7|) |§ has only 
qualitative character as well. 



VII. Conclusion 

Our consideration of the quartet state of the two-dimensional antiferromagnetic Heisenberg 
model was sufficiently consistent. Following to || and [|Hj we have solved exactly the problem 
for four spins and determined the initial ground state and structure of higher excitations. At 
the next step we produced mapping of our quartet state into the Bose system with many degrees 
of freedom. At the low temperatures and level of excitations this Bose system is reduced to the 
system with the nonhermitian polynomial Hamiltonian. We constructed the ground state of 
this Hamiltonian in the mean- field approximation and found a nontrivial solution of the self- 
consistent equations for the normal and abnormal averages. We calculated the basic corrections 
to the mean-field approximation and showed that they are small. We took into consideration 
the most essential higher triplet excitations and calculated their contribution to the ground 
state energy. We also found that these contribution are small. We did not consider the higher 
singlet and quintet excitations but their contribution to the ground state energy definitely small 
because the direct transitions from the initial ground state into these excited states are absent. 
We did not find an unstable excitations in the quartet state. 

AS a result we found that the ground state energy is sufficiently high E /N = — 0.6J 
per spin. We believe that this result has a quantitative character due to the effective small 
parameter of the perturbation theory. It is in our case 1/z where z is the number of spins 
neighbors to the given quartet. In practice, the corrections are small due to smallness of the 
angle integrals over the Brillouin zone. 

Therefore, it is followed from our result, that the quartet state of the two- dimensional 
antiferromagnetic Heisenberg model can not be in competition with the Neel state as it was 
suggested in flOf . Is the quartet short range order in competition with the Neel short range 



order for the doped antiferromagnet? It is an open question at present. 



Acknowledgments 

We would like to thank S.A.Kivelson and S.V.Maleev for the stimulating discussions, A.F.Barabanov 
and L.A.Maksimov for the helpful discussions and sending reprints of their works on the quar- 
tet model. This work was supported partly by the Council on Superconductivity of Rus- 
sian Academy of Sciences, Grant No. 90214 and by the scientific-technical program "High- 



11 



Temperature Superconductivity" as part of the state program "Universities as Center for Fun- 
damental Researches". 



Appendix A: Structure of a single quartet 

Every quartet consists of the four spins and has the Hamiltonian H q : 

H q = J(s a s b ) = (J/2) ((s a + s b ) 2 -si- si) , (A.l) 

where s a = Si + s 3 , s b = s 2 + s 4 . The states of a quartet can be numbered by the four quantum 
numbers: S is the total spin of a quartet S = s a + s b , s a , s b are the spins of a and b subsystems 
and m is the projection of the total spin. The general state is \s a s b Sm > and we have the 
following states: 

\(p >— ]1, 1, 0, >, \tm >— |1, 1, l,m >, \ip >= |0, 0, 0, >, 

\qm >— |1, 1, 2, m >, \am >— |1, 0, 1, m >, \bm >— |0, 1, 1, m > . (A. 2) 

The energies of these states are following: 

E v = —2J, E t = —J, E a = E b = = 0, E q = J. (A. 3) 

It is convenient for our purpose to use linear combinations 

\xm >= (\am > -\bm >)/v / 2, \ym >= (\am > +\bm >)/v / 2 (A.4) 

instead the states \am > and \bm >. In this work we restrict our consideration by the ip -singlet 
and the a, x, y- triplet. The wave functions of this states can be presented in a form: 




|yi>=2 



l.o >= T2 

\y-l>=- 



+ 



+ 
+ 



(A.5) 



where signs ± note the direction of corresponding spins. The state \<p > can be also presented 
as a sum of two possible dimer configurations: 



\<p >= (\dl2 > <8>|d34 > +\dU > <g>|d23 >)/v / 3, 
\dab >= (\a+ > ®\b— > —\a— > ®\b+ >)/V2. 



(A.6) 



This representation shows that our primary state sufficiently close to the RVB-state. 



Appendix B: Structure of selfconsistent equations 



When solving the system of equations (|3.12| ) we transform the two-dimensional over k inte- 



gral into the one-dimensional over the energetic variable e one. This is possible because the 
dependence over k is not explicit but only via 7k : 



1 i 



E/(7k) = / P(e)f(e)de, p(e) = £ 5(e - 7k ), J p(e)de = 1. (B.l) 



k -i k _! 

One can easily get for p(e) 



P(e) = -MO, e = xfi^e, (B.2) 



7T 2 ' 



where K(£') is a complete elliptic integral of the first kind. For the computation of integrals 
over e we have used the Gauss type formula: 



i 

/ 



/0(e)/(e)de = $>/(*)■ (B.3) 

i=i 



Which is precise for only polynomial of the order 2n — 1. The values £j are determined by 
zeros of the orthonormalized polynomials n-th order on the segment (-1,1) with the weight p(e) 
. The coefficients q also can be expressed in terms of these polynomials |T3. Substituting 



the expressions ( |3.6| ) for a^, (3^, 5^ in the terms N p , F p into Eqs. ( |3.12| ) we get the following 
selfconsistent system of the equations 

N p = (3/4 + 4iVi + 4F+) I p + (4iV + F + - 1 + 3^/4) I p+1 - 5 p0 /2, 
F p = -(2N, + 2F 1 ')I P - (8N - 1 - 3Ff /4)/ p+1 , 

F p + = -{2F+ - 1 - 3F+/4)W (B.4) 
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Here the integrals I p are determined by the equations 

f P = / P -E(7k) p J R k 1 (n k + l/2) = ! p = 0,l,2. (B.5) 

k 

As we can see from (B4) the five unknown variables N , Aq, F + , F+, F{~ can be expressed in 
terms of the three integrals I p for p = 0, 1, 2 which are determined by Eqs.(B5). These variables 
I p were found by searching zeroes of the function 

$ = (/o) 2 + (A) 2 + (/ 2 ) 2 - (B.6) 
This method of solving of Eqs.(B4) is more effective than the direct solution of these equations. 
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